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We study the proximity effect in hybrid structures consisting of superconductor and ferromagnetic 
insulator separated by a normal diffusive metal (S/N/FI structures). These stuctures were proposed 
to realize the absolute spin- valve effect. We pay special attention to the gaps in the density of states 
of the normal part. We show that the effect of the ferromagnet is twofold: It not only shifts the 
density of states but also provides suppression of the gap. The mechanism of this suppression is 
remarkably similar to that due to magnetic impurities. Our results are obtained from the solution 
of one-dimensional Usadel equation supplemented with boundary conditions for matrix current at 
both interfaces. 

PACS numbers: 74.45.-(-c, 72.10.-d, 74.78.-w, 75.70.-i 



I. INTRODUCTION 

The research on heterostrusturcs that combine su- 
perconductors and ferromagnets has begun in sixties i 
Still, the FjS structures remain a subject of active ex- 
perimental and theoretical investigation. New devel- 
opments concern Josephson, 7r-j unctions,^ spin valves 
based on giant magnetoresistance efFect^S. triplet super- 
conducting ordering^ Andreev reflection phenomena in 
S — F trilayersS. and proposed detection of electron 
entanglement.® 

Near the F/ S interface electrons are influenced by both 
exchange field h of the ferromagnet and pair potential 
A of the superconductor. Exchange field tends to split 
the density of states so that the energy bands for dif- 
ferent spin directions are shifted in energyi Besides, the 
exchange field at F / S interfaces induces pair breaking, 
suppression of the spectrum gap and even formation of 
a superconducting gapless stateJ. The latter is qualita- 
tively similar to the gapless state in superconductors with 
paramagnetic impurities.^ The exchange field is also ac- 
tive if the ferromagnet is an insulator (FI): Although in 
this case electrons can not penetrate the ferromagnet, 
they pick up the exchange field while reflecting from the 
FI/S interface. ^"-^^ This has been experimentally veri- 
fied with EuO-Al| AI2O3 I Al junctionsi^ 

The physics of F/S or FI/S structures is thus gov- 
erned by two factors: i. electron states are modified by A 
and h, ii. the A and h are modified as a result of the mod- 
ification of electron states by virtue of self-consistency 
equations. A and h correspond to incompatible types of 
ordering that suppress each other and therefore compete 
rather than collaborate. It was suggestedi^ that S/N/FI 
structures can be used to get rid of the second factor. 
The buffer normal metal effectively separates A and h in 
space to prevent their mutual suppression, provided its 
thickness exceeds the superconducting coherence length. 
However, the electrons in the normal do feel both su- 
perconducting and ferromagnetic correlations. Varying 
the conductances of the S/N and N/FI interfaces it is 
possible to tune the strength of these correlations. 



If there are no ferromagnetic correlations, the tradi- 
tional proximity effect in S/N structures takes place. A 
S/N interface couples electrons and holes in the normal 
metal by the coherent process of Andreev reflection^ 
at energies e ~ Aji^ so that Andreev bound states are 
formed.?'*' If the normal metal is connected to the bulk 
superconductor only, there is a (mini) gap in the spec- 
trum of these states. The energy scale of the gap is not 
A. Rather, it is set by the inverse escape time into the 
superconductor, H/te ^ A. This minigap was first pre- 
dicted in Ref. [ijand has been intensively investigated 
A common realistic assumption is that diffusive transport 
takes place in the normal metal. In this case, the su- 
perconducting proximity effect is described by the Usadel 
equation.^" 

The S/N/FI structures can be used to achieve an 
absolute spin-valve effect^ The collaboration of super- 
conducting and ferromagnetic correlations results in a 
spin-split BCS-like DOS in the normal metal part, very 
much like in a BCS superconductor in the presence of 
the spin magnetic field. ^ However, the effective exchange 
field h and proximity gap A characterizing the DOS in 
this case,^« are parametrically different from h and A 
in the ferromagnet and superconductor. In particular, 
the fact that the effective exchange field h affects elec- 
trons in all points of the normal part of a S/N/FI struc- 
ture, does not imply to that the "real" exchange field h 
in the ferromagnet penetrates into N by an appreciable 
distance. Actually it is known that h only penetrates up 
to distances of the Fermi wavelength in the normal part. 
Rather, the effect of h is due to the extended nature of 
the electron wave functions in N , which probe the spin 
dependent potential at the FI /N interface and carry the 
information about this potential throughout the whole 
normal metal region. 

Two such stuctures with normal metal parts connected 
by a tunnel junction constitute the absolute spin valve^ 
The presence of a normal metal is essential to provide a 
physical separation between the sources of superconduct- 
ing and magnetic correlations so that superconductivity 
and magnetism do not compete. The fact that the mag- 
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net is an insulator guarantees the absence of the normal 
electrons at the Fermi level and thus enables the prox- 
imity gap. The spin- valve effect mentioned would not 
be absolute if the ferromagnet is a metal. In this case, 
there would be a finite density of states at any energy 
due to the possible electron escape into the ferromagnet. 
Besides, the use of a insulator does not require nearly 
fully spin polarized ferromagnets (half-metal materials) 
to achieve an absolute spin- valve effect. 

The analysis of Ref. was restricted to the so-called 
"circuit-theory" limit. The variation of Green functions 
along the normal part was disregarded. This is justified if 
the resistance of the normal metal part is much smaller 
than both the resistance of the S/N interface and ef- 
fective spin- mixing resistance characterizing the N/FI 
interface. 

In the present paper, we extend this analysis to arbi- 
trary resistances of the diffusive normal metal part. To 
do so, we analyse the solutions of one-dimensional Us- 
adel equationSS in the normal part. Our goal is to find 
the gaps in the spectrum for both spin directions. The 
equation must be supplemented by boundary conditions 
at both magnetic and superconducting interfaces. 

Microscopic models for interfaces of mesoscopic struc- 
tures have been extensively studied in past years com- 
bining the scattering matrix approach and quasiclassi- 
cal Green's function theoryiSii2SiSi24 In the case of dif- 
fusive conductors the interfaces can be described in a 
compact /transparent way by means of "circuit theory" 
boundary conditionsiS^ In that case the interface is de- 
scribed by a set of conductance parameters given as cer- 
tain specific combinations of the reflection and/or trans- 
mission amplitudes of the scattering matrix associated 
with the interface. The regime of diffusive transport 
considered is distinct from the (quasi)ballistic regime as- 
sumed in many studies of F/S structures liSiii*^^ 

To summarize the resuls shortly, we have shown that 
the effect of magnetic correlations is twofold. Firstly, 
these correlations shift the BCS-like densities of states 
in energy, with shifts being opposite for opposite spin 
directionsiiLH In the first approximation, the absolute 
value of the proximity gap is not affected by the ferro- 
magnetic insulator. Secondly, we have also found that 
the magnetic correlations may suppres the gap. The gap 
completely dissapears at some critical values of the pa- 
rameters. This is qualitatively different from Ref. 13 and 
presents the effect of the resistance of the normal metal. 

The mechanism of the gap suppression appears to 
be surprisingly similar to that due to paramagnetic 
impurities.^ At qualitative level, this has been noticed 
in the context of FI/S structures>i2iii However, for 
S/N/FI structures the analogy becomes closer: the gap 
closing in a certain limit (see Section lyTl is described 
by equations identical to those of Ref. 0. We stress 
that there are no magnetic impurities in our model and 
the quasiparticles are affected by magnetism only when 
they are reflected at the N/ FI boundary. Effective spin- 
flip time thus arises from interplay of magnetic correla- 



tions and diffusive scattering in the normal metal. As we 
have already mentioned, the spectrum gap suppression in 
S/N/FI structures is not accompanied by suppression of 
the pair potential, this is in distinction from the situation 
described m.^-^^-^i 

The structure of the article is as follows. In section IITI 
we introduce the basic equations and define the matrix 
current in the diffusive normal metal. In section UTTI we 
specify the boundary conditions for the Usadel equation 
by imposing matrix current conservation at the interfaces 
of our S/N/FI system. The resulting set of equations al- 
lows us to calculate the total Green's function G at any 
point in the diffusive wire. Analytical solutions can be 
only found in two limiting cases f sections IIVI and lV)l . Fur- 
ther on, we solve numerically the equations for general 
values of the parameters to obtain the general bound- 
ary in parameter space that separates gap and no-gap 
solutions fsection lvH) . We conclude in Section IVlII 



II. MATRIX CURRENT AND USADEL 
EQUATION 

Let us consider a S/N/FI structure with the diffusive 
normal metal in the form of a slab of length L and cross- 
section S. This accounts both for sandwich L <C ^fS and 
wire L » ^/S geometries. The Usadel equation in the 
normal part can be presented as 

Gn ^ {g{x)^G{x)] ^-^^[e f^, G{x)] (1) 



dx 



where Gn = S/L is the conductance associated with 
the diffusive metal, a being its conductivity, Gq = e^/ft, 
S (X 1/SL is the average level spacing in the metal and e 
is the energy of the quasiparticles (electrons and holes) 
with respect to the Fermi energy Ep. G{x) is an isotropic 
quasiclassical Green's function in Keldysh®Nambu®Spin 
space, which is denoted by ( v ). in Eq. (Q, x is the 
coordinate normalized to the length h, x = 0(1) corre- 
sponding to the superconducting (ferromagnetic) inter- 
face. 

It is convenient to define the matrix current^- as 



/ 

Gn 



-Gix) — G{x) 



(2) 



Substituting Eq. ^ into Eq. ^ we present the latter 
as conservation law of the matrix current: 



'-I) = Ileakage{x) Iq - h = I dxlleakage{x) , 

JO 



dx 



(3) 

where "leakage" current Iieakage{x) = —i 
Gq [e fa/^, G(a;)] . Note that the leakage current 
Iieakage{x) does not Contribute to the nonequilibrium 
(physical) charge current, given the fact that its Keldysh 
component is zero. This implies that the physical 
(charge) current is conserved through the system, as 



expected. On the other hand, the diagonal components 
of heakage (x) in Kcldysh spacc (Retarded and Advanced 

components lUi'tagei^) )' which are proportional to the 
energy e, describe decoherence between electrons and 
holes. 

By solving the Usadel equation in the normal metal, 
we obtain the Green's function G{x) that contains 
information about the equilibrium spectral properties 
(G^('^)(x)) and about the nonequilibrium transport 
properties {G^{x)).'^^ In this paper we concentrate on 
spectral properties of the diffusive wire, so from now on 
we restrict ourselves to the retarded block in Keldysh 
space. This is denoted by (^). Note that the retarded 
Green's function G^{x) = G{x) is still a matrix of gen- 
eral structure in Nambu(X)Spin space. 

If there is a single ferromagnetic element in the sys- 
tem, G(x) is diagonal in spin space and can separated 
into two blocks for spin parallel (t) and antiparallel (|) 
to the magnetization of the ferromagnet. For each spin 
component, G{x) can be parametrized in Nambu space 
as G{x) = cos9{x) + sin6'(a;) cos(j){x) fi -I- sin0(a;) 
sm(l){x) f2, fi.2,3 being Pauli matrices. If there is a single 
superconducting reservoir attached to the normal metal, 
the transport properties will not depend on the absolute 
phase of the superconductor, so that (/){x) can be set to 
zero (l){x) = 0. Then G(x) depends on one parameter 
only: G{x) = cos6'(a;) -|-sin6'(a;) fi. The phase (dif- 
ference) (f>{x) maybe important if the normal metal is 
connected to two or more superconducting reservoirs. 

Using this parametrization for G{x), the retarded block 
of Eqs.([3) and Q transforms into a differential equation 
for the angle 9{x): 



—eix)+i ^sm9ix) 



and 



-i G 



de{x) 



N- 



T2 



(4) 



(5) 



where we have introduced Et = hD / L? = Gn S/Gq, 
the Thouless energy associated with the normal metal. 



III. BOUNDARY CONDITIONS 

Circuit theory allows to find the boundary conditions 
at the interfaces of the normal metal in contact with the 
reservoirs simply by imposing matrix current conserva- 
tion at these points. The matrix currents to the reservoirs 
are given by circuit theory expressions. (FIG. 1). We will 
assume that the superconducting reservoir is coupled to 
the normal metal through a tunnel contact. In addition, 
we disregard energy dependence of Green functions in 
the reservoir assuming that the energy scale of interest 
is much smaller than the superconducting energy gap A 
in the reservoir. Under these assumptions, the retarded 
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FIG. 1: Matrix current and boundary conditions. Circuit- 
theory expressions give matrix currents Is, If from the cor- 
responding reservoirs. These currents should match matrix 
currents jo,i from the Usadel equation, defined via deriva- 
tives of Green functions. This fixes the boundary conditions 
for the Usadel equation. 



Green function in the reservoir is just fi. The matrix 
current to the reservoir thus reads 



Gs 



T1,G{X) 



(6) 



To describe the matrix current to ferromagnetic insula- 
tor, we use the results of our previous workSi where we 
obtain 



f -Get, 
Ip = I — - 



M a Tz,G{x) 



(7) 



CT being matrices in spin space, M being the magnetiza- 
tion vector of the ferromagnet. 

The parameter has a dimension of conductance and 
is related to the imaginary part of so-called mixing con- 
ductance G0 =ImG^^. Mixing conductance has been in- 
troduced in Ref. to describe the spin-flip of electrons 
reflected from a ferromagnetic boundary and is, in gen- 
eral, a complex number. For insulating ferromagnet it is 
however purely imaginary. In that case, G^ acts as an 
effective magnetic field. Such spin-dependent scattering 
situation is shown to be important in magnetic insula- 
tor materials^^ and half-metallic magnets.^® Evaluation 
of G^ for a simple model of insulating ferromagnet can 
be found in Ref. 13. 

Using the parametrization in terms of 9{x), we find 



Is = -G5 cos 6*0 T2 



(8) 
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and 

ipi = ±iG^ sin 9i T2, (9) 

where the +(— ) sign corresponds to up (down) direction 
of spin with respect to the magnetization axis of the fer- 
romagnet, 6*0 = 6^(0) and 9i = 61(1). 

Equating Is = lo, Ipi = h gives the boundary condi- 
tions required, 

-5scos0o = ^|o (10) 



±150 sin 6*1 = 1 1 (11) 

where we have intoduced two important dimensionless 
parameters characterizing the stucture: gs — Gs/Gn, 
g<t> = G(f,/GN- As above, ± denotes two spin directions. 

The solutions of the Usadel equation generaUy corre- 
spond to complex 9. The density of states in a given 
point at a given energy is i^(e, a;) = t'oRe[cos(6'(e, x))], i^o 
being the density of states in the absence of proximity 
effect. In this paper, we concentrate on gap solutions 
where v = Q everywhere in the normal metal. 

For this porpuse, it is convenient to define the complex 
angle 9{x) = 7r/2-|-i //(cc). Real ^(x) corresponds to gap 
solution. In terms of this angle, the full system to solve 
reads 

—^li{x)+ecos\i^{x)=0 (12) 




FIG. 2; (a) g^/gs versus /io for various values of gs f Eg 1181 . 
The maximum achieved at /xq gives the maximum g^ at which 
the gap persists at given gs. (b) The maximum g^ versus gs- 
The curve saturates at g,^ ^ 0.663 for grg — > 00. 



-hgssmh/xo = lo 



±30 cosh ^1 = 



(13) 



(14) 



where we introduce the dimensionless energy e — 2 e/ Et- 
The gap solutions exist in a certain region of the three- 
dimensional parameter space {gs, g^,s). To determine 
the boundary of this region is our primary task. 

There are three limits where the solutions can be ob- 
tained analytically. The limit of vanishing resistance of 
the normal metal, gs — — e, can be treated with 
cirquit theory and has been considered in Ref. Below 
we address the zero energy limit (e « 0) and "spin-flip" 
limit (5s-50<l)- 

IV. ZERO ENERGY 

Here we will analyze the gap solutions at the Fermi 
level, that is, at zero energy. In this case the leakage 
current conveniently dissapears, Iieakage (x) = and / = 
constant. Eq. Q can be easily integrated giving Gi — 



exp{I/GN)Go, where Go(i) 
obtain 



G(0(1)). From this we 



arccos ( Tr ( Go Gi 



4- Tr Go Gi 



Gq, Gi 



(15) 



which can be further simplified by using the parametriza- 
tion Go(i) = 0036*0(1) ■?3 ± sin 6*0(1) t'i into the following 
simple expression for the current through the diffusive 
normal metal: 



-i G 



hT2 



(16) 



Taking into account the boundary conditions on both 
interfaces, we obtain the equations for /io.i, 



gs sinh = Mi - Mo = ±50 cosh fii . 



(17) 



We can readily express from these two equations g^/gs 
as a function of /io and gs 



sinh fiQ 



gs cosh {fio + gs sinh ^o ) 



= /(Mo, 9s) 



(18) 



Here we concentrate on the spin-up component. The so- 
lution for spin-down component corresponds to different 
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sign of ^0- In FIG. 2a, we plot g^jgs versus /xq for sev- 
eral values of gs- We see that gci,/gs reaches a maxi- 
mum value {g^/gsT at a certain value /Xq. The position 
and the height of the maximum changes by changing gs, 
^IiOy^oo,{g^/gsr^li^gs^O. 

So for a given 175, the maximum possible value of 
such as there is still a gap in the induced density of states 
of the wire is given by 



V. "SPIN-FLIP" LIMIT 

Now we would like to extend the results of the previous 
section to finite energy e. We do this assuming that the 
Green function does not change much along the normal 
metal, so that \9i — 9o\ <^ do- 

We start again with Eqs. Q), lfTU|l and ((TT|l . Inte- 
grating Eq.l^l over x and using Eas. l|10|) and 1)11(1 gives 



9*4, = f{tJ'0,9s)9s 



(19) 



In FIG. 2b we plot as a function of gs- This curve 
defines the boundary between gap (below) region and no 
gap (above) region in g^ — gs parameter space at zero en- 
ergy e = 0. As expected, magnetic correlations combat 
the proximity effect at the Fermi level and the gap solu- 
tions dissapear upon increasing g^. Even if the coupling 
to the superconductor is infinitely strong, gs — > 00, the 
gap survives only if g^ < 0.663. 

Let us expand g^/gs near its maximum value {gcjilgs)* - 
Defining deviations from this point Ag^-s = {9<j>/9s) ~ 
75)* and fl = fiQ — iiq, the expansion is obviously 



C fL^ 



(20) 



C being a positive constant. Let us note that the density 
of states i^(0)/i'o = Im sinh(/io) cx Imjl. This gives a 
square-root singularity of the density of states near the 
boundary. 



\/^9<t>- 



(21) 



The limit 5s ~* deserves some special consideration. 
It corresponds to the circuit-theory limit where the dif- 
fusive normal metal simply reduces to a "node" with 
spatially independent Green's function G. This gives 
however, a BCS-like inverse square-root singularity in 
the density of states at the boundary given at gs = g,j,, 

vjva = 1/^1 - {9sl94>Y ~ 1/^2^.90-5 at Ag^^s < L 
This result looks difficult to reconcile with the result ob- 
tained in Eq.|(2U. 

The point is that for 55 ^ 1 an extra crossover takes 
place such the density of states changes from inverse 
square root to square root. The crossover occurs close 
to the gap boundary gs — g,jj- From Eq. (fT5|l it is 
clear that for 55 <g; 1, the values of gcf>/gs — 1 occur 
at /xq ^ 1. Close to the boundary, an expansion of Eq. 
(|18|l up to terms ^ exp(— 2/io) is possible. The evalua- 
tion of the maximum of g^ /gs as function of /io allows 

to determine the constant C in Eq. ((20|l . C ~ 9^^- The 

crossover occurs then at Ag^^s — 9s 1- Below the 

2/3 

crossover, at Ag^^s 9s 1 there is a square-root sin- 
gularity v/vq on gg^^'^ Ag^^s that changes to v/vq ~ 
1/ ^2Ag^_s above the crossover, at Ag^_^s ^ 9^"^- The 
maximum density of states is therefore v/vq ~ 9s^^'^- 



±ig,f,smOi + gs cosOq + i s dxsm0{x) — 0. (22) 

"'0 

If we assume that the Green function does not change 
with X, 9{x) — 60, we derive from Eq. H22|) the circuit- 
theory equation 



i (e ± 50) sin 6*0 -I- gs cos 6*0 = 0. 



(23) 



With this, we reproduce the results dicussed in ReffTs! : 
the density of states mimics the one of a BCS supercon- 
ductor in the presence of the spin-splitting magnetic field 



v{e) 



\^±94>\ 



VW^94>) 



9l 



(24) 



The presence of a gap is strictly speaking a non- 
perturbative effect. The perturbation expansion of the 
density of states shown in Eq. H24|l is valid at high en- 
ergies |e ± /i| ^ A. The leading "spin-dependent" cor- 
rection is proportional to {A'^h)/e^ = {9s9<l>)/^^- This 
implies that the leading diagram in perturbation series 
involves four tunneling amplitudes at the S/N inter- 
face and two spin-dependent refiection amplitudes at the 
N/FI interface. 

The equation 123|) is correct in the leading order in 
^t9St94> 1- There can be however a problem if e is 
too close to =F50 since the first coefficient is anomalously 
small in this case. To account for this, we should re-derive 
these equations to the next-to-leading order. 

We present 9{x) in the form that satisfies boundary 
conditions, 

9{x) =9o + {9i~9o)x + 9'^^\x); 9^^Ho) = 0'^^\l) ^ 

(25) 

and evalutate the corrections 9i—9q, 9'^^\x) in the leading 
order. To evaluate 9i — 9q, let us multiply Eq.Q by x 
and integrate by parts the first term dx x 9{x) where 
{9{x) = d'^9{x)/dx'^), obtaining the following expression 

9i ~{9i-9o)+i e [ dx X sm9{x) = 0. (26) 
Jo 

We set 9{x) — 9o under the sign of integral and use Eq. 
((TT|l to obtain 



±9, +2 



sm 9o - 



(27) 
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With the same accuracy, the differential equation for 
O^^^^lx) reads 



£sin((?o) = 0. 



This results in 



0(1) (a;) 



I e—-(x-x ). 



(28) 



(29) 



Finally we substitute Eqs. (j^S)) to Eq. to get the 

following relation: 

i(e ± f/^) sin^o + 3s cos6'o — Css sin^o cos6'o = 0, (30) 

where C, = (^g^^ g^e + /"ij / gs- As we mentioned 

above, C, plays a role only if |e ± g^\ <^ £,50- Under 
these conditions, the energy dependence of C can be dis- 
regarded and C = gl/'igs- 

The relation H3U|) resembles very much one of the most 
important equations in the superconductivity theory that 
was first derived by Abrikosov and Gor'kov^ to describe 
suppresion of superconductivity by magnetic impurities. 
Precise association is achieved by the following change of 
notations: 




C Ss + Ssy 
•V. . ^ 

^ =0 g^=+gs + & 




{e±g^)ET ^ E, 
gsEr -> A, 
CA = ETglli ^ l/r, 

where A,Ts are respectively energy, superconduct- 
ing order parameter and spin-flip time due to magnetic 
impurities^ This is why we refer to the limit under con- 
sideration as to "spin-flip" limit. To remind, there are 
no magnetic impurities in the structure considered, and 
effective spin-flip comes from interplay of diffusion in the 
normal metal and reflection at the N/FI interface. 

Maki has demonstrated that similar equation accounts 
for gapless superconductivity in a variety of circum- 
stances, C being the pair-breaking parameter.'^" We deflne 
9q ~ 7r/2 + i /io, u ~ tanhuo, w = (e ± 5^) /gs, to be 
close to notations of Ref. y^- This gives 

The maximum value of lo with respect to real u, uo* gives 
the energy interval around e — ^g^ where the gap solu- 
tions occur. This value is determined by the condition 
dco/du — 0, which gives 




achieved at m = w*, 

/ , \ 1/2 

u* = (^l- C'/^) . (33) 



FIG. 3; (a) The sketch of the gap domain in e — plane 
at ^ 1. (b) Exact results (squares) for the maximum 
at = 1, follow the expected y/3gs dependence for small gs- 
The lower curve is the one plotted in FIG. 2. 

There are no gap solutions if C > 1- The region where 
these solutions do occur is sketched in FIG. 3a in e — g^ 
coordinates. It looks like a 45° slanted strip, the width 
of the stip in horizontal direction being given by 2uj*gs. 
Near the origin C = 0. In this situation, it is obtained 
from Eas. (|32l I33|l that the width is 2gs and jeitg^l — 
gs. The gap solutions at zero energy e = dissapear at 
9<t> — 9S- The width gradualy reduces with increasing g^ 
due to the increase of (. The strip ends if C = 1; that is, 
at Iff^l = x/3^»ffs (FIG. 3b). 

This demonstrates that even in the limit of g^, gs <^ ^ 
the ferromagentic insulator not only shifts the gap states 
but also reduces and finally suppesses the gap due to 
effective spin-fiip. We show in the next section that the 
same picture is qualitatively valid for arbitrary values of 
these parameters. 

VI. GAP-NO GAP BOUNDARY IN GENERAL 
CASE 

So far we have studied the boundary separating gap 
and no-gap solutions in the parameter space for two lim- 
iting cases that allow for analytic solutions. In this sec- 
tion, we find this boundary for arbitrary values of the 
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FIG. 4: The gap domains (strips) in — e plane calculated 
for different values of = [0.1, 0.2, .., 50, oo]. Their shape 
is similar to the sketch in FIG. 3a. The tip of each strip 
gives the maximum value of at which the gap survives and 
simultaneously the corresponding energy. 




FIG. 5: Side view of the boundary surface in three- 
dimensional (ffs, 50, e) space. The upper curve gives the max- 
imum connecting the tips of the strips plotted in FIG. 4. 
It rises monotonically with gs to reach the asymptotic value 
~ 2.68 at — > oo. The lower curve presents a cross-section 
of the surface with the e = plane. 



parameters. We do this by solving Eqs. (|12I13I14|I nu- 
merically. 

Since we only concern with the boundary, the numer- 
ical procedure is as follows. We fix gs and e. The so- 
lutions of Eq. El with boundary condition (|13f) can be 
parametrized with /j,o. We express ^i, /ii in terms of /iq. 
Then the last boundary condition (|14|l could be solved 
to find /Lto in terms of g^. We do the opposite: We use 
(|14|l to directly express in terms of fiQ and find the 
two extrema of g<;f>(^o)- Those give the endpoints of the 



interval of where the gap solutions exist — elements of 
the boundary. We plot these endpoins at fixed gs versus 
dimensionless energy e to obtain slanted strips similar to 
the one in FIG. 3a. At certain energy, the extrema come 
together indicating the endpoint of the strip. 

In FIG. 4 we show these strips in — e plane for a 
wide range of values of gs- As expected from the previous 
discussion, for gs ^1 the strips extend along the g,p = £ 
line. The sharp tip of each strip gives the critical value 
of g^ at which for a given gs the induced minigap disap- 
pears. For small gs, the height of the tip, ^/'Sgs, is much 
bigger than the width of the strip gg. With increasing 
gs, the shape of the strips changes. They increase both 
in width and height, so that these dimensions become 
of the same order. The strips also become less slanted. 
The shape converges at gs — > cxd (outer curve in FIG. 4). 
In this limit, the maximum g^ that allows for supercon- 
ductivity is w 2.68 and is achieved at e w 4.71. It is 
interesting to note that this energy is higher than the 
maximum value of the minigap without magnetic cor- 
relations {e{g^ = 0) = 7r/2, see FIG. 4). Counterintu- 
itevely, the presence of the magnetic insulator helps the 
gap solutions to persist at higher energy. Albeit the mag- 
netic correlations quickly remove these solutions from the 
Fermi level. 

Each strip is a cross-section of the boundary surface in 
three-dimensional {gtj,,gs,£) space. We present in FIG. 
5 the side view of this surface. The lower curve in this 
figure is the cross-section of the surface with e = plane 
and shows the critical value of at which the gap solu- 
tions dissapear from the Fermi level. The same curve has 
been already presented in FIG. 2. The upper curve gives 
the critical value of g^ at which gap solutions dissapear 
at any energy. It consists of the tips of each strip from 
FIG. 4 as a function of gs- We see that at 55 — > oo this 
curve satures at = 2.68. The asymptotics g^ = \/igs 
derived in the previous section agree with this curve at 
50 < 1 as expected. 



VII. CONCLUSIONS 

We have studied the proximity effect in S/N/FI struc- 
tures with N being a diffusive normal metal. We pay 
special attention to the gap in the density of states and 
find its domain in the parameter space. The convenient 
dimensionless parameters to work with are g^ , gs that 
characterize the intensity of magnetic and superconduct- 
ing correlations respectively, and energy in units of Thou- 
less energy, e. 

We demonstrate that the combined effect of a ferro- 
magnetic insulator and the elastic scattering on the prox- 
imity gap of a diffusive wire is twofold. First, the ferro- 
magnetic insulator provides an effective exchange field h 
that shifts the gap edges in opposite directions for op- 
posite components without reducing the energy interval 
where the gap solutions occur. Second, its effect com- 
bined with suficiently strong elastic scattering in N re- 
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duces this interval and finally suppresses the gap. In the 
limit of small gs,g<f> ("spin-flip" limit) the mechanism of 
this suppression is precisely equivalent to the known one 
from magnetic impurities, with an effective spin-flip rate 
1/ts = g^Ex/i- Qualitatively, this picture remains valid 
at arbitrary parameters. 

If 50 > 0.66 no gap persist at the Fermi level. If g^ > 
2.68 no gap occurs at any energy. Counterintuitively, the 
gap in the presence of magnetic correlations may occur at 
energies higher than in the absence of the ferromagnetic 
insulator. 

The absence or the presence of the gap in the nor- 
mal part of a S/N/FI structure at certain energy can 
be observed by a spin-sensitive tunnel probe. The resis- 
tance of such probe must exceed all interface resistances. 



The possible implementation of the probe depends on its 
concrete geometry. For traditional sandwich geometry, 
it is probably simpler to keep the FI layer rather thin, 
so that the electrons can tunnel through. Covering this 
layer with a conducting ferromagnct makes the tunnel 
probe. For the wire geometry, small tunnel contacts to 
ferromagnets can be made in different points of the nor- 
mal wire. An alternative is the suggestion of Ref. 0: 
the tunneling between two FI/N/S structures. 
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